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Abstract 

We build, using group-theoretic methods, a general framework for 
approaching multi-particle entanglement. As far as entanglement is 
concerned, two states of n spin-1/2 particles are equivalent if they are 
on the same orbit of the group of local rotations ([/(2)"). We give 
a method for finding the number of parameters needed to describe 
inequivalent n spin-1/2 particles states. We also describe how entan- 
glement of states on a given orbit may be characterized by the stability 
group of the action of the group of local rotations on any point on the 
orbit. 
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1 Introduction 



Discovered in 1964 by J. Bell the existence of non-local correlations among 
remote quantum systems is one of the most fascinating quantum phenomena. 
But while for long time these correlations were considered more as a curiosity, 
recently they have found a large range of applications, forming the very bases 
of quantum communication and quantum computation; obviously the interest 
in better understanding these correlations has increased dramatically. 

Traditionally, starting with Bell, the example which has been most stud- 
ied was that of non-local correlations between two remote quantum parti- 
cles. However, it is now clear that the correlations among more than two 
remote particles present novel and highly nontrivial aspects compared to 
two-particles entanglement. (See for example the correlations generated by 
the GHZ state 0.) Nevertheless, very little is known yet about multi-particle 
entanglement. It is the aim of the present paper to take a few steps towards 
understanding the general structure of multi-particle entanglement. 

The key element in our approach is to note that two states which can be 
transformed one into another by local operations (unitary transformations) 
are equivalent as far as their non-local properties are concerned. This leads 
us to investigate the properties of the Hilbert space of n spin-1/2 particles 
under local unitary transformations. 

We find the following picture emerging. 

• Each particular state \E' belongs to an equivalence class comprised of 
by all states which can be obtained from \l/ by acting on it with local 
unitary operators; all states in a class are equivalent as far as non- 
locality is concerned. Obviously, the Hilbert space of states decomposes 
completely into equivalence classes, or "orbits" , under the action of the 
group of local unitary transformations. 

• An arbitrary state \1/ of n spin 1/2 particles is described by 2*^ complex 
parameters. Some of these parameters (or functions of them) spec- 
ify the equivalence class to which \l/ belongs. These parameters (or 
functions) are obviously invariants under local transformations. The 
reminder describe where \I' is situated inside the equivalence class - 
they do change under local transformations. 
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Incidentally, two-particles entanglement is technically so much simpler 
to study than multi-particles entanglement because there is a simple 
way to identify the invariant parameters - the Schmidt decomposition. 
Indeed, let Cj ® Cj, i,j = 1,2 be some arbitrary base vectors in the 
Hilbert space of the two particles then a general state of two particles 
is given by 

^ = Y.a,,e,(g)e,. (1) 

However, by choosing some appropriate base vectors for each particle, 
the double sum in (1) can be reduced to a single sum 

^ = E/5^/*®/- (2) 

i 

The Schmidt coefficients are manifestly invariant under local transfor- 
mations. Indeed, local unitary transformations can only change the 
Schmidt base vectors, but not the Schmidt coefficients. 

^ = J:^^f^ ® ^ ^' = E ® /i- (3) 

i i 

As it is well-known, for multi-particle states in general there exists 
no similar simple decomposition What can one than do? Given 
our above analysis of multi-particle entanglement, it is now clear that 
instead of simply trying to find something which formally resembles 
the Schmidt decomposition, we should try to follow its spirit, not its 
form. That is, to try and find a representation which separates local 
and non-local parameters. 

As an important result, we find that for large n, most of the parameters 
describe non-local properties. This is opposite to the case of small n - 
for two spins, out of the 8 real parameters which describe a generic (un- 
normalized) state, only 1, the unique independent Schmidt coefficient, 
has non-local significance. 

• Finally, we note that in the case of two-particle entanglement some 
of the states are, in some sense, special. Such states are the direct- 
products and the singlet-like states. We show that the special nature 
of these states is determined by their invariance properties. Namely, 



2 



for these special states there are more local actions which leave them 
unchanged than in the case of generic states. For example in the case 
of a singlet 

;^(ei® 62 - 62 <8)ei) (4) 

where Ci and 62 represent spin polarized "up" or "down" say, along the 
z axis, identical rotations of the two spins leave the state unchanged. 

Furthermore, such enhanced invariance properties are in fact common 
for all states in an equivalence class, and thus characterize the class 
itself. To find the "special" equivalence classes, we have therefore to 
study their invariance properties. We argue that these "special" classes 
describe fundamentally different types of entanglement while a generic 
class represents a combination of different types of entanglement. 

Group-theoretically, the situation is the following. The space of states of 
n spin 1/2 particles is the n-fold tensor product (D^" =(D^ ... and the 
group of local transformations is the n-fold product \J (2)" = [/(2) x ... x f/(2), 
(each copy of U(2) acting on a different spin, i.e. on the corresponding 
copy of C^). The equivalence classes are orbits under the action of the local 
transformations group. Hence, the space of orbits is 

C/(2) X ... X C/(2)' 

this is the main mathematical object we are investigating. 

The number of parameters needed to describe the position of ^ on its 
orbit is the dimension of the orbit. Not all orbits have the same dimension. 
As noted above, there are "special" orbits - singular orbits- which have higher 
invariance, i.e. lower dimension. 

The total number of parameters (2" complex parameters = 2"+^ real 
parameters) describing the space of states minus the number of parameters 
describing a generic orbit (the dimension of the orbit), gives the number of 
parameters describing the location of the orbit in the space of orbits, i.e. the 
number of parameters describing the non-local properties of the states. 
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2 The number of parameters needed to de- 
scribe inequivalent states 

In this section we are interested in finding out how many parameters are 
needed to describe the space of orbits of the action of f/(2)" on the space of 
states, i.e. the number of parameters which describes inequivalent states. To 
do this it will be convenient to find the (real) dimension of a general orbit; 
the number of parameters is then found by subtracting this number from 

A lower bound on this number can be obtained by a simple argument 
of counting parameters. Each of the n copies of the local unitary group 
U{2) is described by 4 real parameters. Thus there can be no more than 
4n parameters describing local properties of the states, and hence at least 
2"+i — 477, non-local parameters (i.e., invariants under local transformations)^. 

One can immediately see that for large n almost all parameters have 
non-local significance. 

The above bound is, in general, not satisfied. The reason is that not all 4n 
parameters describing the local transformations lead to independent effects. 
For example, equally changing the phase of all states of any particular spin 
has the same effect as changing the phases of any other. Hence, at least, the 
group of local transformations reduces from f/(2)" to U{1) x SU{2)"' which 
has dimension 3n -|- 1. This leads to a better lower bound on the number of 
non-local parameters of 2"+^ — (3n + 1). 

This is, however, not the end of the story. We will find below that the 
number of parameters describing independent local transformations may be 
fewer (and correspondingly, the number of non-local parameters larger). 

2.1 Dimension of a general orbit 

To find the dimension of a general orbit it is simplest to work infinitesimally. 
Thus, in general, associated to the action of each element of a Lie algebra 
of a Lie group K which acts on a space V there is a vector field: take an 

■^For convenience we always consider non-normalized states, and thus the norm also 
appears as one of the invariant parameters. 
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element T of a basis for the Lie algebra, the action of the group element 
k = exp ieT G -ftT on an element v E V induces an action on functions from 
V to C; and the vector field, Xt, associated to the Lie algebra element T is 
found by differentiating: 

XTf{vf^§-J{e-^v)U=o. (6) 

The linear span of vector fields at the point v associated with the whole 
Lie algebra forms the tangent space to the orbit at the point v and so the 
number of linearly independent vector fields at this point gives the dimension 
of the orbit. 



2.2 A single spin 

The case n ~ 1 helps to illustrate the general formalism. The space of states 
has real dimension four (complex dimension two). It is also clear that the 
action of a unitary operator on a vector cannot change its norm, so that the 
dimension of the space of orbits must be at least one (in fact we will soon see 
that it is precisely one). However, the group U{2) has dimension four so that 
the set of vector fields associated to an arbitrary basis for the Lie algebra 
cannot be linearly independent. 

In the representation of U (2) acting on a convenient Hermitian basis 
for the Lie algebra is 

ay = 

Now take an element 

(8) 

and consider the infinitesimal change under a transformation in the direction 
So that under a group transformation close to the identity. 
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We now write everything in terms of real variables: 

a — ci + idi, (3 — C2 + io?2- 

so that 



(11) 
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(12) 

Thus there is an induced action on a function f{v) — /(ci, di, C2, (^2): 

/(ci, di, C2, (^2) "-^ /(ci - ed2, di + ec2, C2 - edi, d2 + eci). (13) 
Differentiating with respect to e we find: 



de 



6 = 



, d d . d d \ „ 

-d2Tr- + C2Trr - "1^ + cittt J- 
oci odi 0C2 od2 1 



(14) 



We write the vector field associated to this Lie algebra element as 

C2 



,9 9 , (9 d 

-02^ h C2^J- - dl- h Ci — 

oci odi 0C2 od2 



Ux-V where 



-dl 
V ci / 



(15) 

In a similar way we may find the vectors Uy, and Ui associated to trans- 
formations by CTy, and I2: 



( C2 \ 

d2 

V -dl J 



U2: = 



Ui = 



(16) 



/ -dl \ 

Cl 

d2 

V -C2 / 

It is not too difficult to check that only three of these four vectors are linearly 
independent. Indeed 



/ -dl \ 

Cl 

-d2 

\C2 I 



2((ii(i2 + CiC2)Mx + 2(cirf2 — rflC2)My+(c^ + (il— C2— (i2)'2^z~(Ci + C?i + C2 + (i2)Wl = 0. 

(17) 

Thus the dimension of the orbit is three and so there is one parameter (the 
norm) which describes the different orbits. 
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2.3 Two spins 



In a similar way we may analyze the case of two spins. A general vector may 
be written 



(18) 



i,j=l i,j=l 

where {61,62} is a general basis of(D^. In the representation of t/(2)^ on 
we may use the following basis for the eight Lie algebra elements: 

cTj; (g) I2, (8) I2, CFz^h, h ® h,h® CTx, I2 <8) cr^, I2 <8) c^^, I2 ® l2- (19) 

One sees that the clement I2 (E) I2 appears twice, so that in fact there are only 
seven different Lie algebra elements to consider. If we choose the following 
order for the coordinates of the eight dimensional real vector space: 



(Cii, du, C12, di2, C21, (^21, C22, d22), 

then the derivative operator is 

„,ddddddd 



d 



^dcii ' ddii ' dci2 ' ddi2 ' dc2i ' dd2i ' 8022 ' dd22 



(1) ^^(1) 



(20) 



(21) 



y ■<y -2, 



and the vector fields are all of the form ?i. Vg- The vectors {u}^',u, 
u^^\u^\u^P} associated to the Lie algebra elements {a^ I2, CTy 1 
az (8> I2, I2 ® I2, 12 ® cr^, I2 ® cTiy, I2 ® cr^} respectively (the superscript on 
u refers the component in the tensor product, the subscript the Lie algebra 
element) are 



4^^ = 


( — (i21, C21, —^22 


C22, — (ill, cii, — (ii2 
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— Cii, —(ill, ~Ci2, - 
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C12, — (i2i, C21, —d22 






(— (ii2, C12, —dii 


Cll, — C?22, C22, —d2l 




uf = 


(Ci2, (ii2, — Cii, - 


d\i, C22, d22, — C21, - 


-d2ir, 


uf = 


(—(ill, Cii, (ii2, - 


-C12, —(^21; C21, 0?22; " 


-C22f 



(22) 

It may be shown that only six of these vectors are linearly independent for 
general values of the Cij and dij. Thus the dimension of the generic orbit 
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is six and therefore the number of parameters describing the different orbits 
is two. This confirms the well-known result that any state of two spins is 
equivalent, under local rotations, to one of the form 

A^(cos0ei®ei + sin0e2®e2). (23) 
2.4 Three spins 

A computation similar to the one in the above subsections shows that in 
the case of 3 spin 1/2 particles the dimension of a generic orbit is 10, and 
hence the number of real non-local parameters (including the norm) is 6 
(= 23+1 _ 10). 

It is interesting to note that in this case all the 3 x 3 + 1 = 10 parameters 
describing the local transformations U{1) x SU{2Y are actually independent. 

By brute force one can show that any 3 spin 1/2 particle state is equiva- 
lent, up to local transformations to 

cos aei (cos /3ei ® ei + sin /5e2 62) + 
N sin a cos 762 (S> (sin /5ei ® Ci — cos /3e2 62) + 

sin a sin 762 (S> (cos 5ei ® 62 + e^^ sin 6e2 ® ei ) . (24) 

A systematic way of finding the invariants is given in the next section. 

3 Invariants 

For some purposes one might wish to know whether or not two states are 
on the same orbit, i.e. are equivalent. In principle one can take the ideas of 
the previous section further to find invariants of the orbits. For consider any 
function on the space of states. If it is invariant under the action of the group 
then in particular it is invariant under infinitesimal group transformations. 
Thus it must be annihilated by the vector fields associated to the infinitesi- 
mal group transformations. Therefore in order to find a set of infinitesimal 
invariants one has to solve a set of simultaneous partial differential equations; 

^This result was found independently by J. Schlienz [|| 
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the number of such equations is the number of hnearly independent vectors 
associated with the Lie algebra, as in the previous section. 

If we label the Lie algebra elements of local transformations {Tj}, i = 
1...3n+ 1 (corresponding to the local transformations group U{1) x SU{2)^, 
see section 2) then the vector fields X^. are derived as in eq. (||) and an 
invariant function satisfies 

XtJ' = 0, ^ = 1...3n + l, (25) 

a set of 3n + 1 simultaneous linear partial differential equations. The 
method of characteristics allows one to solve the problem in principle, subject 
to being able to perform the integrals which arise. Unfortunately, one can 
easily see that the problem becomes very difficult, even for two spins, for in 
this case one has to solve six simultaneous partial differential equations^. 

It may turn out to be more profitable to realize that one can write down 
a series of polynomial expressions which are manifestly invariant under the 
local actions. We will first show a few examples and then discuss the general 
case. 

3.1 Examples 

In the case of one spin, with general state 

2 

^ = E«^ei, (26) 
1=1 

one can easily see that the expression 

E«^«: (27) 

1=1 

(i.e. the norm of the state) is invariant under local unitary transformations. 

In the case of two spins, with general state \1/ = J2'ij=iC(ij^i ® ^j, the 
norm of the state is invariant and given by a similar expression: 

2 2 

h= J2 (^ij(^hjAii^ni = J2 aij^lj- (28) 

^As we saw in section 2.3 only six of the seven vector fields are linear independent in 
tliis case. 
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There is, however a second, quartic, expression which is functionally inde- 
pendent of /i which is also clearly invariant, since the indices have been 
contracted with the invariant tensor S: 

2 

^2 ^ikO^iim^j'mi^jiki^ih^jji^kki^mmi 
1 

= ^aika*^aj^a*f. = Trace (^(aa'^) ^ . (29) 
In the familiar form of the Schmidt coefficients eq. (^) 

h = N^ 

l2 = N\cos'^(f) + sm^(f)). (30) 

Since we know that in the case of two spins there can only be two invariants, 
any further invariants must be able to be written in terms of Ji and J2. For 
example, consider 

/s = ajfca*^aj„a*fcaprnap„ = Trace Maa^M . (31) 

i,j,k,m,n,p=l 

By noting, for example, that the 2x2 matrix aa^ is hermitian and satisfies 
a quadratic equation (by the Cayley Hamilton theorem), one may show that 

h = l {^hh - I't) . (32) 
In a similar way one may see that all higher order invariants are of the form 

Jjv = Trace , > 3 (33) 

and are expressible in terms of Ji and J2. 

3.2 General case 

A generic state of n spin 1/2 particles can be written as 

2 

^= aiii2...inen ® (g> ... ® ei„. 

il,i2,..i„=l 
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Then a general polynomial expression in the coefficients is 

If the polynomial (|34D has equal numbers of a and a* and all the indexes of 
a are contracted with those of a*, each index being contracted with an index 
located on the same slot (i.e. if Ci-^....k„... are appropriate products of 5's) then 
the polynomial is manifestly invariant. 

For example, in the case of three spins with generic state 

2 

^ = H ctijkei ® Cj (g) Cfe, (35) 

i,j,k=l 

there is one quadratic invariant, the norm, there are the quartic invariants 
(in addition to the square of the norm), 

2 

'Jl — 2^ ^ijkC^ijm'^pqmOipqk 
1 
2 

<^2 — 2^ '^ikjO^imj^pmqOipkq 
1 

2 

•^3 ^kijC^rnij^'rnpqC^kpqi (^^) 

1 

and so on, the different invariants arising by contracting indices in different 
ways. 

Furthermore, one can prove that a// invariant polynomials are constructed 
in this way. The proof of this theorem (not given here) is based on the fact 
that all polynomial functions of k vectors in (D^, invariant under U{2) are 
polynomials in the inner-product of the vectors0. 

A key issue is the following. There are infinitely many polynomial invari- 
ants. We need to be able to construct from these polynomials a complete set 
of functionally independent invariants for arbitrary numbers of spins. For- 
tunately, given any set of polynomials, there is a algorithmic procedure of 
determining the relations between them using the theory of Grobner bases 
[0]. Thus there is a systematic way of constructing sufficient independent 
invariants to classify states; indeed the procedure applies to more general 
groups and representations than U{2) acting on (D^. Firstly determine the 
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number Nj of independent invariants using the ideas of the previous section. 
Then determine the number of independent invariants at lowest order (in the 
case of U{2)^ acting on ®"C^ there was just one, the norm, of the form aa*). 
If this number is less than Nj, construct the invariants at the next order 
and see which of these are functionally independent of each other and the 
ones previously constructed. The procedure continues until Nf are found. 
There is a simple formula, the Molien formula, for the generating function of 
the number of linearly independent invariants at each order which may well 
simplify the task0. 

4 Orbit Types 

As discussed in the introduction, a further important question that the group 
theoretic approach allows one to address is what types of entanglement can 
occur. One can do this by recalling that by definition any group G acts 
transitively on an orbit O and thus an orbit may be written as 

O = G/H (37) 

where H is the stability group of any point on the orbit. Thus the space of 
states of n-spins, (D^ breaks up into orbits each of which is characterized by 
its stability group. Each stability group is a subgroup of U{2)"', so the issue 
is then to find which subgroups occur as stability groups. A generic orbit 
will have a certain stability group, but there are also special cases are where 
an orbit has a larger symmetry group. If we denote by the invariance 
group of the state \E', we will see that states with "maximal" symmetry are 
particularly interesting. By states of "maximal" symmetry, we mean those 
states \l/ for which there are no others which have an invariance group which 
contain if^ as a proper subgroup. 

One systematic way to analyze the space of states, in principle, is to use 
the infinitesimal methods of section 2. Consider the 2-spin case. We found 
that of the eight generators of f/(2)^, only six were linearly independent for 
generic states so that generic orbits have a two dimensional invariance group. 
However there will be some values of the parameters describing the states for 
which the number of linearly independent vectors is smaller than six. Finding 
these points is a problem in linear algebra. Unfortunately the complexity of 
the calculation seems to make it impractical. 
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An alternative approach is to make use of the fact that every stabihty 
group is a subgroup of [/(2)". One can make a hst of subgroups of U (2)" and 
check which subgroups occur as stabihty groups. Goursat's theorem [Q] gives 
a complete characterization of subgroups of any direct product of two groups 
and this enables one, in principle, to produce this list. The complete set of 
subgroups, even of U{2) x U{2) is considerable, once all discrete subgroups are 
taken into account. However, the example below shows that much progress in 
understanding the space of states can be made by considering only continuous 
subgroups in the first instance. 

As an example, consider a (fairly general) three-spin state of the form 

\l/ = aei ® ei ® ei + be2 (S) 62 ® 62 + cei ® ei 62 + de2 ® ei ei. (38) 

In order to find whether this state is invariant under any continuous 
(connected) group, it suffices to check whether it is annihilated by any Lie 
algebra element. As mentioned in section 2, since each copy of U{2) in the 
group [/(2)^ contains a U{1) subgroup corresponding to changing the global 
phase of the state, it suffices to consider SU{2Y x f/(l); thus the phase is 
counted only once. The most general Lie algebra element in this case is 

T = ai(a^)i + a2(0-x)2 + Oi3{a^)s + Pl{(Ty)i + P2icry)2 + P3i(^y)3 

+71(0-2)1 + 72(0-^)2 + 73(0-2)3 + sis, (39) 

where 

(o-^)i = 0-^012® I2; {a^)2 = I2 ® cTx- ® I2 etc. (40) 
and Ig is the identity element 

Ig = I2 ® I2 ® h- (41) 

By direct calculation one can check that if a, b, c and d are all non-zero, then 
the state is not annihilated by any non-zero Lie algebra element so that the 
state is not invariant under any continuous (connected) group. 

The special cases, where the state does have an invariance group, are 
interesting, however: consider first the case a = 0. If 6, c and d are all non- 
zero then we find that the state is annihilated by the Lie algebra element 
with 7i = —72 = 73 = 5 with all other coefficients in T being zero; i.e. the 
state is invariant under U{1). 
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If however a = b = and c and d are non-zero with |c| 7^ \d\, then we find 
that invariance is further enhanced and the state is invariant under U{1)'^. 
If |c| = \d\, the state has yet further symmetry, namely C/(l) x SU{2) and 
one notices that the state is of the form a singlet with respect to particles 1 
and 3 tensor product with a vector for particle 2; we write this as singlet^^ ® 
vector2. The invariance group U{1) x SU{2) arises since a singlet is invariant 
under a (diagonal) SU{2) and the state vector-2 is invariant under U{1). The 
invariance group of the state cannot be increased by choosing special (non- 
zero) values of c and d so a state of the form singleti^ vector2 has maximal 
symmetry. 

If a = 6 = and one of c or c? are also zero, wc find the symmetry is 
also enhanced with respect to the case where c and d are non-zero: in this 
case the symmetry is U{1)^ and such a state also has maximal symmetry in 
the sense that no state has symmetry group of which this is a subset. The 
state is of the form (g> ■u;2 ® "^3 (i-e. it is homogeneous). In the case that 
a = b = c = the generators may be taken to be {(Jz)i + Is, {o'z)2 — Is and 
{c!'z)3 — Is, for example. 

One also finds a similar structure among the states with a — and c — 
or a — and d — 0, namely invariance group of U (1)^ in unless the state is one 
of the special ones with maximal symmetry namely either homogeneous with 
invariance f/(l)^, or of the form singlet® vector invariance SU (2) x t/(l). 

The cases of the sets of states with 6 = or c = have similar structure 
to those with a — Q. The case of d = is different, however. 

If d = and a, b and c are all non-zero, one calculates that the state is 
annihilated by {(yz)i ~ (^'2)2 only; the state is invariant under U{1). If a? = 
and a = but b and c are non-zero, the invariance is enhanced to f/(l)^, in 
general or SU{2) x ^7(1) when \b\ = \c\ in which case the state is of the form 
singleti2 (8) vector^, a state of maximal symmetry. When d — and 6 = 
but a and c are non-zero, the invariance is enhanced to U{lY; the state is 
homogeneous. 

Perhaps the most interesting case is when d = and c = but a and b are 
non-zero, in which case one finds, for all values of a and b, that the state is 
invariant under U{l j^. However although there are a number of states with 
this symmetry, thought of as an abstract group, as described above, the way 
that the group acts on the states is quite different in the case d = c = than 
for example d = a = 0. In the case d = c = 0, the generators are {0-^)1 — {o'z)2 
and {<Jz)2 — (c"z)3; corresponding to correlation between spins 1 and 2 and 
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between 2 and 3. In the case of c? = a = 0, the invariance group arises since 
any vector in (D^ is invariant under U{1) and a generic two particle state is 
also invariant under f/(l). 

Amongst those states with d = c = 0, there are some which larger sym- 
metry groups than f/(l)^. If a = or 6 = 0, then the invariance group is 
U{1)^] the state is homogeneous. However the case a = b, while not having 
further continuous symmetry is picked out by the fact that only this state has 
a discrete symmetry of Z2 corresponding to the operation of simultaneously 
flipping all spins. This is the famous GHZ state. 



5 Conclusion 

In this paper we have started to build a general framework for understand- 
ing multi-particles entanglement. Obviously we have taken just a few steps 
here, and there are far more questions still open than answered. For exam- 
ple, it is known that in case of two-particle entanglement, to get a deeper 
understanding of entanglement one needs to take into account not only local 
unitary transformations but also measurements and classical communication 
between the two observers situated near the two particles. Also one has to 
consider actions taken on a large number of copies of the state and not only 
on a single copy as considered here. Nevertheless, it is clear that any "mea- 
sure of entanglement" for multi-particles must be a function of the invariants 
described here. 
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